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1. Introduction 



The study of multi-component system of nonlinear partial differential equations (PDEs) 
is current interest in the theory of nonlinear waves. Of particular concern are the multi- 
component generalizations of the nonlinear Schrodinger (NLS) equation because of their 
wide applicability in real physical systems such as nonlinear optics, nonlinear water waves 
and plasma physics [1-3]. The integrable two-component NLS system has been introduced 
for the first time by Manakov to describe the propagation of polarized electric field in an 
optical fiber and explored in detail by means of the inverse scattering transform (1ST) 
method [4]. After this remarkable work, various variants of integrable multi-component 
NLS systems have been proposed and analyzed by the exact method of solution such 
as the 1ST and Hirota's direct method. In this paper, we consider the following multi- 
component system of nonlinear PDEs which is a hybrid of the coupled NLS equation and 
coupled derivative NLS equation 

= 0, 0- = l,2,...,n), (LI) 

X 

where qj — qj{x,t) (j = 1,2, ...,n) are complex- valued functions of x and t, /i and 7 are 
real constants, n is an arbitrary positive integer and subscripts x and t appended to qj 
denote partial differentiations. The integrabihty of the above system has been estabhshed 
by constructing the Lax pair and an infinite number of conservation laws [5] . For the two 
special cases 7 = and /i = reduced from the system (1.1) which correspond to the 
multi-component NLS and multi-component derivative NLS equations, respectively, their 
integrability has already been verified in [6, 7]. In the context of plasma physics, the two- 
component system with fj, — is a model equation for the propagation of polarized Alfven 
waves. The single bright sohton solution to this system has been obtained by means of the 
1ST [8] . The two-component system with // 7^ and 7 7^ has been derived as a model 
for describing the propagation of ultra-short pulses in birefringent optical fibers, together 
with its soliton solutions [9]. Quite recently, we obtained the general bright multisoliton 
solution (i.e., the A'"-soliton solution with N being an arbitrary positive integer which 
vanishes at infinity) of the two-component system by using a direct method [10]. It is 
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important to remark that the constant ji must be nonnegative to support smooth bright 
sohtons. The negative case is worth studying as well which will be treated in a separate 
issue. 

The purpose of this paper is to extend the results obtained in [10] to the general 
n-component system. Specifically, we present the bright A^-soliton solution of the n- 
component system (1.1) in the form of compact determinantal expressions. Although the 
construction of the solution can be done following the similar procedure to that developed 
in [10] for the two-component system, we provide a novel proof of the solution using the 
expansion formulas for the bordered determinant. 

This paper is organized as follows. In section 2, we first transform the system (1.1) 
to a gauge equivalent system and then recast it to a system of bilinear equations by 
introducing appropriate dependent variable transformations. Subsequently, the bright N- 
soliton solution to the bilinear equations is presented. It has a simple structure expressed 
in terms of certain determinants. In section 3, we introduce some notations associated 
with the bright A^-soliton solution and then prove several key formulas for determinants. 
In section 4, we perform the proof of the bright A^-soliton solution using an elementary 
theory of determinants in which Jacobi's identity and related formulas will play a central 
role. In section 5, we provide an alternative expression of the bright A'^-soliton solution 
and compare it with the corresponding solution presented in section 2. We then find a 
simple relation between two types of solutions by employing the properties of the Cauchy 
matrix. This result leads to an alternative proof of the solution in a very simple way. In 
section 6, we discuss a (2+l)-dimensional nonlocal modified NLS equation arising from 
the continuum limit n — > oo of the system (1.1). Specifically, we demonstrate that its 
bright iV-soliton solution can be generated simply from that of the n-component system. 
Section 7 is devoted to concluding remarks where we will comment on existing literatures 
about the bright A^-soliton solutions of the multi-component NLS equation (equation 
(1.1) with 7 = 0) and show that the solutions obtained in this paper include these known 
solutions as special cases. 
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2. Bilinecirization and bright A^^-soliton solution 



2.1. Bilinearization 

We first apply the gauge transformations 



Qj — Uj exp 



17 
' 2 



\uk\'^dx 



k=l 



(j = l,2,...,n), 



(2.1) 



to the system (1.1) subjected to the the boundary conditions qj — >■ 0,Uj {j — 
1,2, ...,n) as \x\ — >■ 00, where Uj — Uj{x,t) {j — l,2,...,n) are complex-valued functions 
of X and t. Then, we obtain the system of nonhnear PDEs for Uj 



iuj,t + Uj,^^ + n \ukfj Uj + i7 uluk,xj uj = 0, (j = 1, 2, n), (2.2) 

where the asterisk appended to Uk denotes complex conjugate. This notation will be used 
frequently hereafter. The second step in our analysis is given by the following proposition: 



Proposition 2.1. By means of the dependent variable transformations 

% = y, (i = l,2,...,n), 



(2.3) 



the system of nonlinear PDEs (2.2) can he decoupled into the following system of bilinear 
equations for f and Qj 



{iDt + Dl)gj-f^O, (i = l,2,...,n), 

n 
k=l 

n . n 

Dlf ■ f* ^ + 'iJ2^^9k - gl- 



k=l 



(2.4) 
(2.5) 

(2.6) 



k=l 



Here, f — f{x,t) and Qj — gj{x,t) {j — 1,2, ...,n) are complex-valued functions of x and 
t and the bilinear operators and Dt are defined by 



(2.7) 



where m and n are nonnegative integers. 



Proof. Substituting (2.3) into (2.2) and rewriting the resultant equation in terms of the 
bihnear operators, equations (2.2) can be rewritten as 

-i(iA^, • / + Dig, • /) + i-rDlf • / + /x/ l^.r + i7 E alD.gk • / j = 0, 

(j = l,2,...,n). (2.8) 

Insert the identity 

fDlf ■ f = fDlf . r - 2f,DJ . r + f{DJ . f%, (2.9) 
into the second term on the left-hand side of (2.8). Then, equations (2.8) become 

^ n n 

-{iD,g, . f + Dig, ■ f) + ^ [f{-Dlf • /* + //^ " (^-/ " + ^^J^dldk,.} 



k=l k=l 



n 

+U[2DJ ■ r -i^Y.\9k?}] =0' (i = l,2,...,n). (2.10) 

k=l 

As easily confirmed by a direct calculation, the left-hand side of (2.10) becomes zero by 
virtue of equations (2.4)-(2. 6). □ 



It now follows from (2.3) and (2.5) that 

-TZ^I^^I =^^^T' ^2.11) 



fc=i ■' 

which, substituted into (2.1), yields the solution of the system (1.1) in the form 

= (j = l,2,...,n). (2.12) 

Note that for the n-component NLS equation (the system (1.1) with 7 = 0), the solution 
(2.12) simplifies to q, — g,/ f. Indeed, if 7 = 0, then the bilinear equation (2.5) reduces 
to D,ji.f ■ /* = 0. Thus, the ratio /*// turns out to be an arbitrary function of t which 
can be set to 1 under appropriate boundary condition. 

2.2. Bright N-soliton solution 
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We now state the main result in this paper: 



Theorem 2.1. The bright N-soliton solution of the system of bilinear equations (2.4)- 
(2.6) is given by the determinants f and gj {j — 1, 2, ...,n) where 



f = 



A I 
-I B 



9j 



A I 
-I B 0^ 
-a* 



, (i = l,2,...,n). 



(2.13) 



Here, A, B and I are N x N matrices and z, aj and are N -component row vectors 

defined below and the symbol T denotes the transpose: 

1 z^zl 



A — {ajk)\<j,k<N, djk 



= exp(pjX + ipjt), 



R _ N A _ (/^ + nPk)Cjk _ 

-D - [Ojk)l<j,k<N, Ojk - . , Cjk - OisjOisk^ 

Pj ^ P'' s=l 

I = {^jk)i<j,k<N, : N X N unit matrix, 
z ^ {zi,Z2, ...,zn), a.j ^ {aji,aj2, ...,ajN), = (0, 0, 0). 



(2.14a) 

(2.146) 

(2.14c) 
(2.14ci) 



The above bright A'^-sohton solution is characterized by (n + 1)A^ complex parameters 
Pj{j = 1,2, ...,n) and asj{s = 1, 2, ..,n; j = 1, 2, .... N). The former parameters determine 
the amplitude and velocity of the solitons whereas the latter ones determine the polar- 
izations and the envelope phases of the solitons. The conditions pj + 7^ for all j and 
k and pj ^ Pk for j ^ k may be imposed to assure the regularity of the solution. In the 
special case of n = 2, (2.13) and (2.14) reproduce the bright A^-soliton solution presented 
in [10]. The proof of Theorem 2.1 will be given in section 4. 

To simplify the proof of theorem 2.1, the following observation is useful: 



Proposition 2.2. // we introduce the gauge transformations 



f ^f, Qj = exp 



i{^x+Ui\ t 



7 



7 



2// 



~9o, (i = l,2,...,n), 



X — X I ij ^ t — i ^ 

7 



(2.15a) 
(2.156) 



then the bilinear equations (2.4)- (2.6) recast to 

{iDi+Dlrg,-f = 0, U = l,2,...,n), 
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k=l 



k=l 



respectively. 



Proof. The proof can be done by a straightforward calculation. 



(2.16) 
(2.17) 

(2.18) 
□ 



Thus, the form of equations (2.4) and (2.5) is unchanged whereas equation (2.6) be- 
comes a simplified equation with n = 0. Consequently, the proof of the A'^-soliton solution 
may be performed for the corresponding solution with /i = 0. Hence, in the analysis 
developed in the following sections, we put /i — without loss of generality. 

3. Notation and some basic formulas for determinants 

In this section, we first introduce the notation for matrices and then provide some basic 
formulas for determinants. 

3.1. Notation 

We define the following matrices associated with the A'"-sohton solution (2.13) with (2.14): 

A I 
-I B 



D 



(3.1) 



L'(a*;b) 




A I z^' 
D{a.*-z) = I -/ B 0^ 
a* 



D{z*;z) 



A I z^^ 
-I B 0^ 
z* 

(3.2) 

Note the position of the vectors a*, b, z and z* in the above expressions. The matrices 
which include more than two vectors will be introduced as well. For example, 

/A I z^\ /A I z^ z'^^ 

-I B 0^ 0^ 



L'(a*,z*;b,z) 



-/ B b^ 0^ 
a* 
\z* / 



, D(a*,z*;z,z') 



a* 

Vz* y 



(3.3) 



3.2. Formulas for determinants 

Let A — {ajk)i<j,k<M be an M x M matrix with M being an arbitrary positive integer, 
Ajk be the cofactor of the element ajk and a, b, a.j and hj {j — 1,2,..., n) be M-component 
row vectors. The following well-known formulas are used frequently in our analysis [11]: 

M 



d_ 
dx 



\A\ 



j,k=l 



A 

b z 



M 



\A\z- ^ Ajkajbk, 
j,k=i 

|A(ai,a2;bi,b2)||A| = |>l(ai; bi)P(a2; bs)] - |>l(ai; b2)P(a2; bi)|. 



(3.4) 



(3.5) 
(3.6) 



The formula (3.4) is the differentiation rule of the determinant and (3.5) is the expansion 
formula for a bordered determinant with respect to the last row and last column. The 
formula (3.6) is Jacobi's identity. 

The following two formulas may not be well-known but are very important in our 
analysis. In particular, formula (3.7) gives rise to the expansion formulas for the bordered 
determinant (see (3.9) and (3.10) below): 

|A(ai;bi)| ■■■ |A(ai;b„) 



|^(ai, ...,a„;bi, ...,b„)||A 



n-l 



|^(an;bi 



1^(3^) b,, 



(n > 2), (3.7) 



l^l + I^C-e)"" Yl |^(a,i,...,a,^;b,„...,b, 

m=l l<si<...<Sm<n 

1^1 + 5^^^ Yl |^(a.„...,a,„;b,,,...,b,^)| 



ml 

m=l Sl,...,Sm=l 



(3.8) 



Here, e is an arbitrary parameter, the notation bja^ on the left-hand side of (3.8) repre- 
sents an M X M matrix whose (j, A;) element is given by Psj<^sk and n' = min(n, M). The 
formula (3.7) is a variant of the Sylvester theorem in the theory of determinants. 

Proof of (3.7). The proof proceeds by a mathematical induction. For n = 2, (3.7) 
reduces to Jacobi's identity (3.6). Assume that formula (3.7) is true for n — 1 (n > 3). 
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Let L be the left-hand side of (3.7). Recall that the determinant changes its sign if any 
two rows (or columns) are interchanged. Applying this rule repeatedly to L when cin and 
b„ are shifted in front of ai and bi, respectively, 

L = |i(ai,...,a„_i;bi,...,b„_i)Pp-S 

where A — A(a^; b„) is an (M + 1) x (M + 1) matrix which is assumed to be nonsingular, 
i.e., \A\ ^ 0. In view of the inductive hypothesis, L can be written as 

|i(ai;bi)| ••• |i(ai;b„_i)| 



\A 



n-l 



M|n-2 



|A(an_i;bi)| ••• |A(an_i;b„- 
It follows from Jacobi's identity (3.6) that 

|i(a,-;bfc)||^| = |^(a„,aj;b„,bfc)||^| = |^(a„; b„)||^(a,-; bfe)| - |^(a,-; b„)||^(a„; bfe)|. 

which, substituted into L, recasts L into the form 
1 



|^(^n; b^ 



ln-2 



|A(a„;b„)||A(a,;bfc)| - |A(a,; b„) P(a„; b^) 



l<j,k<n-l 



Referring to the property of the bordered determinant, the above expression simplifies to 



|^(a„;b„)|"- 



|^(a„;b„)||A(ai;bi 



|A(a„;b„)||>l(ai;b„_i 



|A(ai;b, 



|A(a„;b„)||A(a„_i;bi)| ••• |A(a„; b„)||A(a„_i; b„_i)| |A(a„_i;b^ 
|^(a„;bi)| ••• |A(a„;b„_i)| 1 

Extract the factor |/l(a„; b„)| from the jth row (j = 1, 2, ...,n — 1) and then multiply the 

last column by the same factor. Then, the resultant expression is seen to be equal to the 

right-hand side of (3.7). □ 

Proof of (3.8) For n = 1, it follows by using the property of the bordered determinant 
that 

Oil + e/3iiQ;ii • • • OiM + e/5iiQ;iM 



A + ebf ai 



Oil ■ ■ • aiM 



-an 



A hi 

-ai e^-^ 
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Repeated use of the above modification yields 



A + e^ bf a. 



A bf 
-ai 



n 

■ 

; 0^ ■■. 0^ 

-a„ 



Expanding the determinant in powers of e ^, it is found that 

n n 

yl + 6j]bfa, =|^| + e"^(-ir E |^(a.„ a,^; b,„ b^Jle'^'^-). 

s=l m=l l<si<...<Sm<n 

The above expression coincides with (3.8) for n < M since n' = min(n, M) — n. For 
M + 1 < n, on the other hand, the determinant |A(asj, a^^; b^^, b^^)! becomes zero 
identically for M + 1 < m < n, as confirmed by the Laplace expansion of the determinant 
with respect to the last m rows, for example. This implies that the summation with respect 
to m is truncated at m = M which is in accordance with (3.8) since n' — min(n, M) — 
M. The second line of (3.8) follows from the facts that any permutation of the indices 
{si, S2, Sm} does not alter the value of the determinant |v4(asi, a^^; b^^, hs^)\ and 
the total number of the permutation is m!, and if the determinant includes at least two 
same rows (or columns), then it becomes zero identically. □ 

Suppose that \A\ ^ 0. Expanding the determinant on the right-hand side of (3.7) with 
respect to the first column and using (3.7) with n replaced by n — 1, we then obtain an 
expansion formula 

1 " 

|A(ai, a„; bi, b„)| = — ^^{-\y'^\A{s.j, bi)||A(ai, a^.i, a^+i, a„; ba, b„)|. 

I I j=\ 

(3.9) 

Similarly, the expansion with respect to the first row gives 

1 " 

|A(ai, a^; bi, bn)| = E(~^)^~^l^(^i' bj)ll^(a2, a„; bi, bj_i, b^+i, bn)|. 

I I j=\ 

(3.10) 

The above two formulas will be used effectively to prove the bilinear equations (2.5) and 
(2.6). 
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4. Proof of the bright A'^-soliton solution 



In this section, we show that the bright A'"-sohton solution (2.13) with (2.14) satisfies the 
system of bihnear equations (2.4)-(2.6). The proof can be performed for // = 0, as noted 
at the end of section 2. We first prove some formulas associated with the determinants / 
and gj {j — 1,2, ...,n) and then proceed to the proof. 

4.1. Formulas 

In terms of the notation introduced in section 3.1 (see (3.1) and (3.2)), / and gj are 
written in the form 

f^\D\, gj^-\D{3i*;z)\, (j = 1, 2, n). (4.1) 

The differentiation rules of / and gj with respect to t and x are given by the following 
formulas: 

Lemma 4.1. 



ft^-'-{\D{z*;z,)\-\D{zl;z)\}, (4.2) 



2 

f,^~\D{z*;z)\, (4.3) 



= {|i^(z*;z,)| + \D{zl;z)\}, (4.4) 

gj,t = - |L'(a*; z,) | + ^ |L>(a*, z*; z, z,) | , (4.5) 

gj,,^-\D{a*;z,)\, (4.6) 

9j,xx = -|-D(a*;z^^)| + ^\D{a*,z*;z^,z)\. (4.7) 

Here, zt, z^ and z^x are N -component row vectors given by Zt = {iplzi,iplz2, ■■.,ip%ZN), 
Zx = {piZi,P2Z2, ...,PnZn) and z^^ = {p\zi,plz2, ...,p%zn), respectively. 

Proof. We prove (4.2). Let D = {djk)i<j,k<2N be a 2A'" x 2A'" matrix and Dj^ be the 
cofactor of the element djk- It follows by applying the formula (3.4) to the determinant 
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/ given by (2.13) that 



. N 



j,k=l 
. N 

2 

j,k=i 



where in passing to the second hne, use has been made of the relations pj^j = Zj^xi Pk4 — 
zl ^. Referring to the formula (3.5) with z = and taking into accoimt the notation (3.2), 
the above expression reduces to the right-hand side of (4.2). A key feature in the proof 
is that the factor {pj + pl)~^ in the element ajk has been canceled after differentiation 
with respect to t. Using formulas (3.4) and (3.5) as well as some basic properties of 
determinants, formulas (4.3)- (4.7) can be proved in the same way. □ 

The complex conjugate expressions of /, fx and gj can be expressed as follows: 

Lemma 4.2. 



D 



A I 
-I B-i-fC I ' 



f:--l\D{z*;z)\, 
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|D(z*;a,-) 



Proof. We prove (4.8). It follows from (2.14a) and (2.14b) that A* 



(4.8) 

(4.9) 
(4.10) 
A^ and B* = 



B^ — ijC^ where C is an x matrix with elements cjk defined by (2.146). These 
relations lead to the expression of /* 

r = 

Since = |A| for any square matrix A, the above expresion reduces to to the right-hand 
side of (4.8) after multiplying the jth row and kth column {j,k = 1,2, ...,N) by a facter 
— 1. Differentiating (4.8) with respect to x and applying formula (3.5) with z = to the 



A* I 




A^ I 


-I B* 




-I {B-i-fCf 
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resulting expression, formula (4.9) follows immediately. The proof of formula (4.10) can 
be done in the same way. □ 
The following formulas will be used in the proof of (2.5) and (2.6): 

Lemma 4.3. 



|D| = |L>| + -|L>(z*;z)|, 
\D{hl;i)\ = \D{a.l;z)l 
|L>(a*;b,) = -|L'(b*;afe)| - ^|L'(b^ z*; a^, z)|, 

|D(a*;z.)| = |L>(b*;z) + ^|L>(b^z*;z,z)|. 



(4.11) 
(4.12) 
(4.13) 

(4.14) 



|D(z*;z)|=2i7^|D(b*;a,)|, 



fc=i 



|5(z*;z)| = -2i7j]p(a^;bfe)|, 



k=l 



where z and hk are N -component row vectors given respectively by z 
and hk = (a^iP^, Q!fe2P2> OikNP*N)- 



(4.15) 



(4.16) 

{zi/pi,Z2/p2, ■■■,Zn/Pn) 



Proof. First, we prove (4.11). A direct calculation using the elements of B and C given 
by (2.14b) reveals that 

Ojk - ^ICjk = Ojk- 

Pk 

The determinant \D\ from (4.8) is now modified to the form 



\D\ 



A 



-I 



I 

B 



where the last line of the above expression follows immediately from the property of the 
determinant. The definition (2.14a) of ajk now gives 



ZiZ, 



Pk _ 

Pj ^Pj 
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In view of the property of the bordered determinant, \D\ is modified in the form 



\D\ 



A I 

-I B 
^ 1 



which is seen to coincide with (4.11) by applying formula (3.5). The proof of (4.12)-(4.14) 
can be done in the same way. Hence, we omit the proof. 

Let us now proceed to the proof of (4.15). To this end, it is to be noted that the 

determinant / can be rewritten in the form 



/ 



A I 
-I B 



where A and B are N x N matrices defined by 

A = {ajk)i<j,k<N, 5,jk 

B = {bjk)i<j,k<N, bjk = 



Pi+Pk ^ 



Indeed, the above expression of / is derived from / from (4.1) by extracting the factors 
Zj and zl from the jth row and kth column, respectively ior j,k — 1,2, N and then 
multiplying the (iV+j)th row and (iV + A;)th column by the factors z* and Zk, respectively 
for j, k = 1,2, N. Using the formulas (3.4) and (3.5) gives an alternative expression of 

fx 

n 

/. = -i7$^|i?(b:;a,)|. 



k=l 



The formula (4.15) follows by comparing this expression with (4.3). The formula (4.16) 

comes from the complex conjugate expression of (4.15). □ 

4.2. Proof of (8.4) 

The proof of (2.4) proceeds following the same procedure as that of the same equation 
forn = 2 [10]. Let Pi be 

P^ = {iD, + Dl)g,.f. (4.17) 
Substituting (4.1)-(4.7) into (4.17), Pi becomes 



Pi = -\D{a*,z*;z,z,)\\D\ + |D(a*;z)||D(z*;z, 



-|^(a*;z,) 



\D(z*:z) 
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-{i|L'(a*;z,)| + |i^(a*;z,,)|}. (4.18) 

Referring to Jacobi's identity (3.6) and the fundamental formula a; |Z) (a; bi)|+/3|D(a; b2) | = 
|-D(a; Q;bi + /3b2)| {q:,/3 G C), Pi simplifies to Pi = —\D{a*; izt + Zxx)\- Since izt + z^x = 
by (2.14a), the last column of the determinant consists only of zero elements, implying 
that Pi = 0. □ 

4.3. Proof of (2.5) 

The equation to be proved is P2 = 0, where 

n 

P2 = ^./-r-^$]kr. (4.19) 

fc=l 

Substituting (4.1), (4.3) and (4.8)-(4.10) into (4.19), P2 becomes 

1 1 . n 

P2 = --\D\\D{z*;z)\ + -\D\\D{z*;z)\ + ^ J] |D(a- z)| |D(z*; a^l- (4.20) 

k=l 

Further simplication is possible with use of (4.11), (4.15) and (4.16) with (4.13), giving 
rise to 

n 

fc=i 

Applying Jacobi's identity (3.6) to the middle term and replacing |£)(b^;z)| by the right- 
hand side of (4.12) in the resultant expression, P2 reduces to 

n 
k=l 

It now follows from (3.8) that 

|P>(a^;z)| = |P>(a:;z)| + ^^ ^ |P>(a^ a^^, a^^; z, a,„ a,J|, (4.23a) 

m=l fei,...,fcm=l 

|5(z*;a,)| = |i^(z*;a,)| + ^^ |i^(z*, a^,, a^^; a,, a,„ a,^)|, (4.236) 

m=l ki,...,km=^ 
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where n" — min(n — l^N — 1). Referring to the expansion formulas (3.9) and (3.10), one 
has 



m 

+p|-^J](-l)'|D(a:^;z)||L>(a*,a*^,...,a*^_^,a*^^^,...,a*^;a,„...,a,^)|, (4.24a) 

1=1 

P(z*,afc,,...,a^^;afe,afei,...,afc^)| = \D\-^\D{z*;a.k)\\D{a.l^, ...,a.ljak,, ...,^kj\ 

m 

+ l^r'$^(-l)'l^(z*;afej||i^(a^^,...,a^^;afc,afc,,...,afe,_,,afe,^,,...,afc^)|. (4.246) 
z=i 

By introducing (4.23) into (4.22) and then using (4.24), P2 takes the form 



' ' m=l 1=1 



X [-|^K;z)IP(z*;afe)||i^(a;^,a^^,...,a;_^,a^^^^,...,a^^;afe,,...,afc^)| 

+l^(afc;z)||r>(z*;afcJ||D(a^^,...,a;^^;afe,afe,,...,afe,_,,afc,^,,...,afc^)| . (4.25) 

Interchange the indices k and A;; in the first term and then shift the row vector a^^ in front 
of a^j^j and the column vector a^ in front of a^^ , respectively. This leads to the following 
relation 

\^\^k-i^kx-i ■■■i^ki-i-i^ki+i-i ■■■>afe^;afci, ...,afc^j| 

~^ \^\^ki-i^kv ■■■■> ^ki-i^^ki+i^ ■■■■> afc„; afei, afc,_^, afc, afe,_^^, a.k^)\ 
~ l^i^ki^ ^km'^ a/c, a/c^, afc,_j, a.ki^i, afc^)|. 

Note that the value of the determinant is not altered since the total signature caused by 
the above manipulation is (— l)2('-i) — 1. Thus, the first term on the right-hand side of 
(4.25) coincides with the second term and cosequently, P2 — 0. □ 

4.3. Proof of (2.6) 
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Instead of proving (2.6) directly, we differentiate (2.5) by x and add tfie resultant expres- 
sion to (2.6) and then prove the equation P3 = 0, where 

n 

P3 = - -'^Yl (4.26) 

k=l 

This reduces the total amount of calculations considerably and the proof becomes trans- 
parent. It now follows from (4.1), (4.3), (4.4), (4.6) and (4.8)-(4.10) that 

1 1 ■ n 

P3 = -2 {|^(z*;z.)| + |i?(z:;z)|} |D|--|D(z*;z)||D(z*;z)|+^ J] |D(a*;z,)||D(z*;a,)|. 

fe=i 

(4.27) 

Differention of (4.15) with respect to x gives 

n 

|D(z*;z,)| + |i?(z:;z)| = -i7 |I?(b*, z*; a^, z)|. (4.28) 

fc=i 

Inserting (4.15) and (4.28) into (4.27), P3 can be put into the form 

n 

^3 = ^E{l^ll^(b^'^*;^fc'^)l + l^(^*;^)ll^(^fe;t.fc)| + |D(a*;z,)||D(z*;afe)|}. (4.29) 

k=l 

Note from (4.11), (4.13), (4.14) and Jacobi's identity (3.6) that 

|D||D(b*,z*;afe,z)| + |D(z*;z)|p(a^;bfe)| 

= -|L>(z*;a,)|{|L>(b^;z) + i|L>(b^z*;z,z)|} 

--|Z^(z*;a,)||5(a*;z,)|. (4.30) 
After substituting (4.30) into (4.29), P3 becomes 

n 

^3 = f + P(a^;Zx)IP(z*;afe)|}. (4.31) 

k=l 

This expression reduces to (4.22) if one replaces by z. Hence, the proof of the relation 
P3 = completely parallels that of P2 = with P2 from (4.22). □ 
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5. Alternative expression of the bright A^-soliton solution 

Here, we persent an alternative expression of the bright A^-soliton solution in terms of the 
determinants with smaller sizes when compared with those given by (2.13). Explicitly, 
we write it as a theorem: 

Theorem 5.1. The determinants f and g'j {j — 1,2, ...,n) given below satisfy the system 
of bilinear equations (2.4)-(2.6): 

f^\A + B'\, g] = 

where A! and B' are N x N matrices and y and a'j are N -component row vectors defined 
below: 

1 V V* 

^ = («jfe)i<j,fe<iv, a^-fc = ^ ^, \ , Vj = exp(5^x + \(^t), (5.2a) 

^ Hj + Hh 

rtl /-Ll \ il^l^~ '^'~1'ik)^jk I II* /r nL\ 

yi yfc ,5=1 

y = (yi, y2, -, Vn), a-'j = {a'ji, aj2, ■■; «jjv)- (5-2c) 

Here, qj (j = 1,2,...,A^) and a'^j {s = 1,2, ...,n;j = 1,2, ...,A^) are complex parameters 
characterizing the solution. 

Proof. The proof of the solution can be performed in the same way as that of (2.13) 
with (2.14). Indeed, the proof of (2.4), (2.5) and (2.6) reduce respectively to the relations 
(4.18), (4.22) and (4.31) in which the matrix D may be replaced simply by the matrix 
A' + B'. □ 

Let us show that the determinants / and gj from (2.13) are closely related to the 
determinants /' and g'j given by (5.1). The following lemma is useful for this purpose: 

Lemma 5.1. The determinants f and gj given by (2.13) can be rewritten in the form 
f^\I + AB\, g^^ 



A' + B' 



a. 



y 





T 



(i = l,2,...,n). 



(5.1) 



I + AB 
-a* 



(i = l,2,...,n). 



(5.3) 
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Proof. Multiplying / from (2.13) by a factor 



B -I 
I O 



and performing the operation 



of matrix multiplication, the first expression of (5.3) follows immediately. Similarly, the 

B -I 0^ 

second expression is obtained if one multiplies gj by a factor 



/ O 0^ 
1 



. Indeed, 



A I 




B 


-I 




I + AB -A 


-I B 




I 







O I 



\I + AB\, 



A 







/ 
B 



0^ 





B -I 0^ 
1 



I + AB 
O 



T 



-A z 
I 0^ 




I + AB tF 



□ 



Since the value of each factor multiplied from the right is 1, (5.3) follows. 

We now establish the following theorem: 

Theorem 5.2. Under the parameterization qj — —p* {j = l,2,...,iV) and a'^^ ~ 
—asj/{2c*) (s = l,2,...,n;j — 1,2, N), the determinants f,f',gj and g'j satisfy the 
relations 

f = c\A\f, 
9j = c\A\gj, (j = l,2,...,n), 

where 



(5.4) 
(5.5) 



N 



Ilm=lipi-Pmy 



1=1 



(Z = 1,2, ...,7V). (5.6) 



The parameters pj (j = 1,2, ...,N) are assumed to satisfy the conditions Pi + Pm ^ fof 
all I and m and pi ^ Pm for I ^ m. 

Proof. Let A be a Cauchy matrix of the form A — {\ p+p* ^ ■ Then, A — {zj5jk)A{z*^5jk). 
Since \A\ ^ by virtue of the well-known formula for \A\ and the conditions imposed 
on the parameters pj (j = 1,2, ...,A^), the inverse of A exists, implying that A~^ exists 
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as well. Actually, it reads A ^ — {z*5jk) ^{zjSjk) ^. Using the explicit expression of 



, I.e., 



P*i+Pk 



[11], the inverse matrix A ^ can be written in the form 

1 



2c* Ck 



Applying the basic properties of determinants to / and Qj from (5.3) gives 



f=\A\\A-' + Bl 



9j = 1^1 



A- 



B A-^7^ 



(j = l,2,...,n). 



(5.7) 

(5.8) 
(5.9) 



The jth element of the column vector A ^z-^ is 



1=1 

2c* A 1 



2r* ^ X^rn=\kPl 



(P/ - Pr) 



(5.10) 



By Euler's formula, the sum in the last hne turns out to be 1 and hence (^4 ^z^) ^ = 2c* / z* 
Introducing this relation into g^^ 



9j 



A- 



B 



-a: 







(i = l,2,...,n). 



(5.11) 



where z = (2c^/2;^, 2c'^/ 2c*^/ z*^) is an iV-component row vector. 

The next step is to modify the determinants /' and g'y By means of the parametriza- 
tion Qj — — p*, Uj from (5.2a) is related to Zj from (2.14a) by the relation yj — Zj~^. 
Similary, the relation c'^f. — Cjk/{4:CjCk) follows from (2.14b) and (5.2b) and the parame- 
terization a'gj — —asj/{2c*). Substitution of these relations into /' gives 



/' 



+ 



2q,+ql 

n£i(4crQ) 

c-^\A-^ + B\, 



(/^ - hqpc'jk 
ij + il 

2c*Ck 1 
p*+Pk Z*Zk 



+ 



{H + i'ypk)cjk 
P*i+Pk 



(5.12) 
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where in passing to the second hne, the factor 1 / (2c*) has been extracted from the jth row 
{j — 1, 2, .., N) and the factor — l/(2cfe) from the kth column {k — 1, 2, .., N), respectively. 
The similar procedure applied to g'j leads to the expression 



9i = c 



-1 



A' 



-a: 



B 





(j = l,2,...,n). 



(5.13) 



The relation (5.4) follows from (5.8) and (5.12) whereas the relation (5.5) follows from 
(5.11) and (5.13). □ 

Thus, we have obtained the two different expressions for the bright A'"-soliton solution 
of the system of nonhnear PDEs (2.2). Explicitly, they read Uj — Qj/f — g'j/f {j — 
1,2, ...,n). 

The following proposition provides an alternative proof of theorem 5.1: 

Proposition 5.1. If f and gj given respectively by (5.4) and (5.5) satisfy the system of 
bilinear equations (2.4)-(2.6), then f and g'j satisfy the same system of equations, and 
vice versa. 

Proof. Substituting (5.4) and (5.5) into (2.4) and using the definition of the bilinear 
operators, 

cMl'W- . f + Dig'. . f) + c'(Dl\A\ . \A\)g'.f' = 0. (5.14) 

The Cauchy type determinant \A\ can be modified, after extracting the factor Zj from jth 
row (j = 1,2, ...,N) and the factor z^. from kth column {k = 1,2, ...,N), respectively, in 
the form 

N 



\A\ = fi{ziz:) f 
1=1 ^ 



1 1 



— exp 



N 



N 



^{pi + p*i)x + i ^{pf - pf)t 



.1=1 1=1 
Differentiation of |A| with respect to x gives 

N 



1 1 



(5.15) 



N 



(5.16) 



=1 
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It immediately follows from (5.16) that 



Di\A\-\A\^2{\A\\Al,-\A\i)^Q. (5.17) 

It is seen from (5.14), (5.17) and the relation c|y4| ^ that /' and g'- satisfy the bilinear 
equation (2.4). The remaining part of the proposition can be proved in the same way 
if one uses (5.17) and the reality of \A\^ i.e., \A\* = \A\ which is a consequence of the 
Hermitian nature of the matrix A. The proof of the converse proposition proceeds in the 
same way if one uses the relation • \A\~^ = in place of (5.17). □ 

6. A continuum model 

The n-component system (1.1) yields a continuum model when one takes a limit n — )■ oo. 
It represents a (2+l)-dimensional nonlocal modified NLS equation of the form 

iqt + qxx + fJ^i^J \q\'^dy^ q + i-f (^J \q\^dyq^ =0, q = q{x,y,t). (6.1) 

Recall that when 7 = 0, this equation reduces to a (2+l)-dimensional nonlocal NLS 
equation proposed by Zakharov [12]. The exact method of solution for equation (6.1) 
can be developed following the same procedure as that for the system of nonhnear PDEs 
(1.1). Hence, we summarize the main results. 
First, application of the gauge transformation 



q — u exp 



17 '•^ '•^ 



2 



u{x, y, t)\'^dxdy 



00 d —00 



u^u{x,y,t), (6.2) 



to the system (6.1) subjected to the the boundary conditions g — >■ 0, — >■ |x| — > 00 
transforms (6.1) to a nonlocal nonlinear PDE for u 

iuf + + fJ^ (^J \u\'^dy^ u + i'j u*Uxdy^u — 0. (6.3) 

The proposition below is an analog of proposition 2.1: 

Proposition 6.1 By means of the dependent variable transformation 

u^^, (6.4) 
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equation (6.3) can he decoupled into the following system of bilinear equations for f — 
f{x,t) and g ^ g{x,y,t) 

OA + Dl)g • / = 0, (6.5) 
DJ ■r='if \9\'dy, (6.6) 

poo • f'CO 

Dlf-f* = fx \g\'dy + ^ / D,g- g*dy. (6.7) 



Proof. The proof proceeds exactly as that of proposition 2.1. Formally, one may simply 
replace the sum Ylk=i integral dy. □ 



It follows from (6.2), (6.4) and (6.6) that 



Q 



gf* 
f 



2 ' 



(6.8) 



which is just a continuum limit of (2.13). 

The following theorem can be derived from a continuum limit of the bright iV-soliton 
solution given by theorem 2.1 and theorem 5.1: 

Theorem 6.1. The system of bilinear equations (6. 5)- (6. 7) admits the following two 
different expressions f, g and /', g' for the bright N-soliton solution : 



f 



A I 
-I B 



+ g' 



A I 
-I B 
-a* 

A' + B' 



0^ 




(6.9) 



(6.10) 



Here, A and B are N xN matrices given respectively by (2.14a) and (2.14b) with Cjk being 
replaced by J^^aj{y)al{y)dy, A' and B' are N x N matrices given respectively by (5.2a) 
and (5.2b) with c'^y, being replaced by J^^aj{y)a'i,*{y)dy and a = a(y) = {ai,a2, ■■■,ckn) 
and a' = a'(y) — (cki; ck2) '^jv) ^'^s N-component row vectors where aj and a'j {j — 
1,2, ...,N) are continuous functions of y. 
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Proof. The proof can be done in the same way as that of theorem 2.1 and theorem 5.1. 
□ 



Theorem 6.2. Under the parameterization qj — —p* and a'j — 
1,2, N), the determinants f,f',g and g' satisfy the relations 



«./(2c*) (i 



/ = c\A\f, 



(6.11) 



9 = c\A\g', 



(6.12) 



where c is defined by (5.6) and the parameters pj {j — 1,2, ...,N) are specified such that 
Pi + P^ ^ foT all I and m and pi ^ p^ for I ^ m. 



Proposition 6.2. /// and g given by (6.9) satisfy the system of bilinear equations (6.5)- 
(6.7), then f and g' given by (6.11) and (6.12) satisfy the same system of equations, and 
vice versa. 



7. Concluding remcirks 

In this paper, we have obtained the two different expressions for the bright A^-sohton 
solution of an n-component modified NLS equation and found a simple relationship 
between them. We also have presented the bright A^-sohton solution of a continuum 
model arising from the system as the number n of the dependent variables tends to in- 
finity. These solutions include, as special cases, existing solutions for a multi-component 
NLS equation. Actually, when 7 = 0, the system of nonlinear PDEs (1.1) reduces to 
an n-component NLS equation. It admits the bright A^-soliton solution of the form 
^3 — 93I f — d'jl f — ..,n). This fact follows from (2.12) and the relations f* — f 
and /'* = /'. See (2.13) and (5.1). The solution (2.13) with 7 = has been obtained in 
[13] using a direct method whereas the solution (5.3) with 7 = has been constructed 
by means of the 1ST [14]. An alternative expression (5.1) of the solution with 7 = has 



Proof. The proof parallels theorem 5.2. 



□ 



Proof. The proof is completely parallel to that of proposition 5.1. 



□ 
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been derived by employing a method of algebraic geometry [15]. For a continuum model 
(6.1) with 7 = 0, the solution takes the form q — g/f. The bright A'"-sohton solution (6.9) 
with 7 = has been found in [16] by a direct method. 

In a future work, we will investigate the various features of the multi-component 
bright solitons. In particular, wc will be concerned with the effect of the parameters ^ 
and 7 which characterize the different types of nonlinearities on the interaction process 
of solitons. 
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